CERN-TH/98-360 



Duality symmetry of Reggeon interactions in multicolour QCD 

L.N. Lipatov * 
Petersburg Nuclear Physics Institute, 
Gatchina, 188350, St. Petersburg, Russia 



Abstract 

The duality symmetry of the Hamiltonian and integrals of motion for Reggeon interactions 
in multicolour QCD is formulated as an integral equation for the wave function of compound 
states of n reggeized gluons. In particular the Odderon problem in QCD is reduced to the 
solution of the one-dimensional Schrodinger equation. The Odderon Hamiltonian is written 
in a normal form, which gives a possibility to express it as a function of its integrals of 
motion. 
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1 Introduction 



The hadron scattering amplitude at high energies y/s in the leading logarithmic approxi- 
mation (LLA) of the perturbation theory is obtained by calculating and summing all con- 
tributions (g 2 ln(s)) n , where g is the coupling constant. In this approximation the gluon is 
reggeized and the BFKL Pomeron is a compound state of two reggeized gluons [1] . Next-to- 
leading corrections to the BFKL equation were also calculated [2], which gives a possibility 
to find its region of applicability. In particular the Mobius invariance of the equation valid 
in LLA [1] turns out to be violated after taking into account next-to-leading terms. 

The asymptotic behaviour oc s JO of scattering amplitudes is governed by the j'-plane 
singularities of the t-channel partial waves fj(t). The position of these singularities ujq = 
jo — 1 for the Feynman diagrams with n reggeized gluons in the t-channel is proportional to 
eigenvalues of a Schrodinger-like equation [3]. For the multicolour QCD N c — ■> oo the colour 
structure and the coordinate dependence of the eigenf unctions are factorized [4]. 

The wave function f m p~2, fti; Po) of the colourless compound state O m ~(p~o) 
depends on the two-dimensional impact parameters ~p~i, ~p~2, p^ of the reggeized gluons. It 
belongs to the basic series of unitary representations of the Mobius group transformations 

apk + b 

Pk^ —j, 1 

cpk + a 

where p k = x k + Wk, p\ = x k — Wk and a, b, c, d are arbitrary complex parameters [1]. For 
this series the conformal weights 

m = 1/2 + iv + n/2, m = l/2 + iu- n/2 (2) 

are expressed in terms of the anomalous dimension 7 = 1 + 2iv and the integer conformal 
spin n of the composite operators O m ~ (p~o). They are related to the eigenvalues 

M 2 f m ,^ = m(m - l)/ m ~ , M* 2 f m ~ = m{m - l)f m ~ (3) 
of the Casimir operators M 2 and M* 2 : 

(n \ 2 

J2 M£) = E 2 K M s = - E Prsdrd s , M* 2 = (M 2 )*. (4) 
k=l ) r<s r<s 

Here M% are the Mobius group generators 

Ml=p k d k) M k - = d k: Mt = -p\d k (5) 

and d k = d/(dp k ). 

The wave function f ~ satisfies the Schrodinger equation [4]: 

^m,rnfm,m ^ f m,m' (^0 

Its eigenvalue E m ~ is proportional to the position ~ = j — 1 of a j-plane singularity of 
the t-channel partial wave: 
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governing the n-Reggeon asymptotic contribution to the total cross-section a tot ~ s m ' m . 

In the particular case of the Odderon, being a compound state of three reggeized gluons 
with the charge parity C = — 1 and the signature Pj = —1, the eigenvalue oj^~ is related 
to the high-energy behaviour of the difference of the total cross-sections a pp and cr p p for 
interactions of particles p and antiparticles p with a target: 



(3) 

u pp u pp ° 



(8) 

The Hamiltonian H in the multicolour QCD has the property of the holomorphic sepa- 
rability [4]: 



H=-(h + h*), [h,h*] = 0, 
where the holomorphic and anti-holomorphic Hamiltonians 



(9) 



h — ^2 h ktk+1 , h* — ^2 K,k+i 

k=l k=l 

are expressed in terms of the BFKL operator [4] : 



(10) 



h k ,k+i = log(pfc) + log(p fc+ i) + — log(p ktk+1 )p k + — !— \og(p kM+l )p k+l + 27. (11) 

Pk Pk+i 

Here p k ,k+i = p k - p k+1 , p k = id/(dp k ),p* k = id/(dp%), and 7 = is the Euler 

constant. 

Owing to the holomorphic separability of h, the wave function / m ^(/0i, pi, fti] Po) has 
the property of the holomorphic factorization [4]: 

/m,m(Pl, pi, ■■; K; Po) = °r,l fm(P^ P2, Pn, Po) f~(P*l,P*2, Pn5 Po) > ( 12 ) 

r,l 

where r and / enumerate degenerate solutions of the Schrodinger equations in the holomor- 
phic and anti-holomorphic sub-spaces: 

e m fm = h f m , e~ /~ = h /~ , E m ~ = € m + e~ . (13) 

Similarly to the case of two-dimensional conformal field theories, the coefficients c r j are 
fixed by the single-valuedness condition for the function f m ~(p~i, pi, p^; po) in the two- 
dimensional ~/?-space. 

There are two different normalization conditions for the wave function [4], [5]: 



J r=l 



d 2 p r 



II Pr,r+1 / 



r=l 



r n 

=/n 



<i 2 p r 



r=l 



(14) 



compatible with the hermicity properties of H. Indeed, the transposed Hamiltonian U is 
related with h by two different similarity transformations [5]: 



h* = n Pr h n Pr 1 = n ^,r+i ^ n ^+1. 



(15) 



r=l 



r=l 



r=l 



r=l 
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Therefore h commutes 



[h,A] = (16) 

with the differential operator [4] 

A = Pl2P23---PnlPlP2---Pn • (17) 

Furthermore [5], there is a family {g r }of mutually commuting integrals of motion: 

[q r ,q,]=0, [q r ,h]=0. (18) 

They are given as 

Qr y ' Pii«2 Pi2'i3"-Piril Pil Pi2'"Pir' (19) 

ii<»2<---<V 

In particular q n is equal to A and q 2 is proportional to M 2 . 

The generating function for these integrals of motion coincides with the transfer matrix 
T for the XXX model [5]: 

n 

T(u)=tr(L 1 (u)L 2 (u)...L n (u)) = J2u n ~ r qr, (20) 



r=0 



where the L-operators are 



«->-( m -^T — a*) — (i 1 >*- (21) 

The transfer matrix is the trace of the monodromy matrix t(u) [6]: 

T{u) = tr (t(u)), t(u) = L l (u)L 2 {u)...L n (u) . (22) 
It can be verified that t(u) satisfies the Yang-Baxter equation [5], [6]: 

qiu) q(v) iZliv -u) = - u) t%{v) tii{u) , (23) 

where l(w) is the L-operator for the well-known Heisenberg spin model: 

ll[ll(w)=wS a s }8tl + i8 s a }8 a s ?. (24) 
The commutativity of T(u) and T(v) 

[T(u),T(v)] = (25) 

is a consequence of the Yang-Baxter equation. 
If one will parametrize t(u) in the form 

t(u)=( j ^+{^ .,^,1 (26) 
this equation is reduced to the following Lorentz-covariant relations for the currents j^u): 
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\jn(u)Ju(v)] = [j^v),j u (u)] = 



% e 



fiupa 



(f(u)f(v)-f(v)f(u)). 



(27) 



2(u - v) 

Here e Mi , pcr is the antisymmetric tensor (ei23o = 1) in the four- dimensional Minkovski space 
and the metric tensor g^ u has the signature (1,-1,— 1,-1). This form follows from the 
invariance of the Yang-Baxter equations under the Lorentz transformations. 

The generators for the spacial rotations coincide with that of the Mobius transformations 



k=i 



Ml = p k d k) Ml= l -{l-pl)d k) Ml= l -{l + pl)d k . 
The commutation relations for the Lorentz algebra are given below: 

M s , M*j = ie stu M u , [ikP, iV*j = it stu N u , [N s , N'] = ie stu M u , 

where iV* are the Lorentz boost generators. 

The commutativity of the transfer matrix T(u) with the local hamiltonian h [5], [7] 



(28) 
(29) 

(30) 



[T(u),h]=0 



is a consequence of the relation: 



[L k (u) L k+1 (u),h k>k+1 ] 



-i (L k (u) 



J k+1 



;«)) 



(31) 



(32) 



for the pair Hamiltonian h ktk+ i. In turn, this relation follows from the Mobius invariance of 
h ktk+ i and the identity: 



h 



k,k+l: 



( *\ 2 

[Mk,k+1 J j N kjk+ i 



AN, 



k,k+l 



(33) 



where 



Mi 



k+l 



M k)k+1 = M k + M k+1 , N k>k+1 = M k 

are the Lorentz group generators for the two gluon state. 

Because the pair hamiltonian h kjk+ i depends only on the Casimir operator (M fe)fe+1 ) , it 
is diagonal 



(34) 
n 2 



h k , k +i \m kM i) = (ip(m k)k+ i) + - m kjk+ i) - 2 ^(1)) \m k ,k+i) 
in the conformal weight representation: 

->\2 



[M kik+1 J \m kik+1 ) = m Kk+l (m Kk+l - 1) \m kjk+1 ) 



(35) 



(36) 



Using the commutation relations between M kik+ i and N kik+ i and taking into account that 

/ >\2 ' — > 

(M k \ = 0, one can verify that the operator N kjk+ i has non- vanishing matrix elements only 
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between the states {mk^+i) and {rrik^+i ±1). It means that the above identity for hk t k+i is 
a consequence of the known recurrence relations for the -^-functions: 

ip(m) =ip(m-l) + l/(m-l) , - m) = i[>(2 - m) + l/(m - 1). (37) 
The pair Hamiltonian hk t k+i can be expressed in terms of the small-w asymptotics 

L k:k+1 (u) = P k ,k+i(^ + iuh ktk+1 + ...) (38) 

of the fundamental L-operator L kjk+ i(u) acting on functions f(p k ,p k+ i) [6]. Here Pk,k+i is 
defined by the relation: 

Pk,k+i f(p k ,Pk+i) = f(pk+i,Pk)- (39) 
The operator L k>k+1 satisfies the linear equation [6] 

L k (u) L k+1 (v) L k)k+1 (u - v) = L k ,k+i{u - v) L k+1 (v) L k (u). (40) 
This equation can be solved in a way similar to that for h kik+ i 



Lk,k+i(u) ~ Pk,k+i» 



T(m k ,k+i + m)r(l - m kt k+i + iu) T(l - iu) 



^ r(m fc)fc+ i - iu)T(l - m k:k+1 - iu) T(l + iu) ' 
where the integral operator m^k+i is defined by the relation 



(41) 



->\2 



m k ,k+i (m k:k+1 - 1) = [M k + M k+1 j (42) 

and the proportionality constant, being a periodic function of m ktk+ i with an unit period, is 
fixed from the triangle relation 

L 13 (u) L 23 {v) L 12 (u -v) = L 12 (u - v) L 23 (v) L 13 (u) . (43) 

To find a representation of the Yang-Baxter commutation relations, the algebraic Bethe 
anzatz is used [6]. To begin with, in the above parametrization of the monodromy matrix t(u) 
in terms of the currents jn(u), one should construct the pseudovacuum state |0) satisfying 
the equations 

j + (u) |0> = 0. (44) 

However, these equations have a non-trivial solution only if the above L-operators are regu- 
larized as 

Li(u)=( U + P^~ i5 P* ) (45) 

by introducing an infinitesimally small conformal weight <5 — > for reggeized gluons (another 
possibility is to use the dual space corresponding to 5 = — 1 [7]). For this regularization the 
pseudovacuum state is 
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n 

\S) = II Pf (46) 

k=l 

It is also an eigenstate of the transfer matrix: 

T(u) \S) = 2j (u) \6) = ((« - i6) n + (u + 1 5) n ) \5). (47) 

Furthermore, all excited states are obtained by applying the product of the operators 
j-(v) to the pseudovacuum state: 

\viv 2 ...v k ) = j-(yx) j-(v 2 )...j-(v k ) \5). (48) 
They are eigenfunctions of the transfer matrix T{u) with the eigenvalues: 

u — v r — i , . r ,„ JL u — v r + i 



T( u ) = (u+ isr n + (« - ity n , (49) 

r=1 u-v r ^ u-v r 

providing that the spectral parameters v i, v 2 , v k are solutions of the set of Bethe equations: 



V s +10 I " V s - V r + I 



for s — 1, 2.../c. 

Due to above relations the function 



k 



Q(u) = U(u-v r ) (51) 



r=l 



satisfies the Baxter equation [6,7]: 



T(u) Q{u) = {u- i5) n Q(u + i) + (u + i5) n Q{u - i) , (52) 

where T{u) is an eigenvalue of the transfer matrix. Its corresponding eigenfunctions can be 
expressed in terms of the solution Q^ k \u) of this equation as follows [7] 

\viv 2 ...v k ) = Q^iu,) QW{u 2 ) ...Q {k \u n ^) \5) , (53) 
where the integral operators u r are zeros of the current j-(u): 

n—l 

j_(u)=cH(u-u r ). (54) 

r=l 

Eigenvalues e of the holomorphic Hamiltonian h also can be expressed in terms of Q(u) [7]. 
Up to now the Baxter equation was solved only for the case of the BFKL Pomeron (n = 2). 
This is the reason why we use below another approach, based on the diagonalization of the 
transfer matrix. 



6 



2 Duality of Reggeon interactions at large N c 



The differential operators q r and the Hamiltonian h are invariant under the cyclic permu- 
tation of gluon indices i — > i + 1 (i — 1, 2...n), corresponding to the Bose symmetry of the 
Reggeon wave function at N c — > oo. It is remarkable that above operators are invariant also 
under the more general canonical transformation: 

Pi-ij -> pi -> p iji+ i , (55) 

combined with reversing the order of the operator multiplication. 

This invariance is obvious for the Hamiltonian h if we write it in the form: 



h = h p + hp , (56) 



where 



K = it f ln (Pk) + Pk,k+x ln(pfe) p J +A + 7 ) , 

k=l V A=±l / 

h P = it ( ln (Pfe,fc+l) + ^ S ^fc+(l+A)/2 Hpk,k+l) Pfe+(l+A)/2 + 7 ) ■ (57) 
fe=l V A=±l / 

Here 7 = — is the Euler constant. 

The invariance of the transfer matrix can be verified using two equivalent representations 
for q r : 

r I h+l \ r I H+l \ 

& = yi n I £ Pfc-i.* p<« ) = z! n p^a+i £ p* ) ■ ( 5s ) 

ii<i2<---V i=l \fc=i(+l / ii<ii<...i r 1=1 \ k=ii+l J 

Note, that the supersymmetry corresponds to an analogous generalization of transla- 
tions to super-translations. Furthermore, the Kramers- Wannier duality in the Ising model 
and the popular electromagnetic duality E <-> H can be considered as similar canonical 
transformations [8]. 

The above duality symmetry is realized as an unitary transformation only for the van- 
ishing total momentum: 

n 

7? = £^ = °- (59) 

r=l 

In this case one can parametrize gluon momenta in terms of momentum transfers k r as 
follows 

p r = k r — k r +i , (60) 
which gives a possibility to present the symmetry transformation in a simpler form: 

k r — > p r — > k r +i , r — 1, 2. ..n . (61) 

Because the operators g r compose a complete set of invariants of the transformation, the 
Hamiltonian h should be their function 

h = h(q 2 ,q3,-,q n ), (62) 
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fixed by the property of its locality. Furthermore, a common eigenfunction of q r (r = 2, n) 
is simultaneously a solution of the Schrodinger equation, which means, that the duality 
symmetry gives an explanation of the integrability of the Reggeon model at N c — > oo. 

To formulate the duality as an integral equation we work in the two-dimensional impact 
parameter space ~ft, initially without taking into account the property of the holomorphic 
factorization of the Green functions. The wave function Vmm of the composite state with 
7? = can be written in terms of the eigenfunction / m ~ of a commuting set of the operators 
q k and q* k for k = 1, 2. ..n as follows 

^m,m(Pi2>P23>->A2) = J ^ / m ~ (pt, p* n \ pt>) . (63) 

It is a highest-weight component of the Mobius group representation with the quantum 
numbers m = 1/2 + iu + n/2 and m = 1/2 + iv — n/2 related with eigenvalues of the Casimir 
operators 

(n \ 2 / n \ 2 

E^ ^m,m = m ( m ~ !)VVm ' E ^fc ^m,S ~ %S ■ ( 64 ) 

k=l I \k=l ) 

The other components of this highest-weight representation can be obtained by applying 
to ip m y n the Mobius group generators: 

/ n \ ri / n \ r 2 

Note that, in accordance with the relations 

m* = 1 — m , m* = 1 — m (66) 
the conjugate function is transformed as /i_ mj i_^: 

(/m,m(Pl>->^;^)) ~/l-m,l-m(Pl>->^;^)- i 67 ) 

Moreover, because of the reality of the Mobius group, the complex-conjugated representa- 
tions are linearly dependent: 

(/ m ,m(pt>--^;Po)) ~/ ^PO' (POO') 2m ~ 2 (Pi5o') 2m ~ 2 / m ,m(pt ) --,P^;PCK)- (68) 

By considering the limit p - >■ oo of this equation, we obtain for ip m> ~ the new represen- 
tation 

^ m ,m(P^> Pil> Aa) ~ /l-m,l-m(Pl> ^> -> °°) ■ ( 6 9) 

Because of the relations 



, m (Pl2>P23,-,P^l)) ~^l-m,l- m (Pi2»P23»-"»Art)» ( 7 °) 

the functions ib~ and ■?/> ~ have the same conformal spin n = m — m. 

T m,m r m,rn r 

Taking into account the hermicity properties of the total Hamiltonian: 
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h + = n \PKk+ii 2 H\{\p Kk+1 \ 2 = n \p k \ 2 H]j \ Pk \- 2 , (71) 

k=l k=l k=l k=l 

the solution ijji of the complex-conjugated Schrodinger equation for ~j? = can be ex- 
pressed in terms of ip~ m as follows : 



^~ m (Pl2,P23,-,Au) = II |P*,fc+il fe m (Pl2,Pi3,-,^l)J • (72) 



fc=l 



If one performs the Fourier transformation of to the momentum space 

n-1 J2 ' 



j Pfc-1, 

with substituting the arguments 



^^^■■■^)=/n-^n^^± m (^S,---^) (73) 



pt -> Pi>K , (74) 
the new expression ^ m ~ (p~u, P23, Pni) will have the same properties as the initial function 
m(Pi2) P2I' •••Pni) under rotations and dilatations. 

Moreover, in accordance with the above duality symmetry it satisfies the same set of 
equations as ip mt ~ and therefore these two functions can be chosen to be proportional: 

^m,m(Pl2>^3>->Aa) = C m,m *m,m(Pl2> P23' ■ (75) 

The proportionality constant c m ~ is determined from the condition that the norm of the 
function 



, n-l ^2 



= /'ff% ±1 <^ (76) 
J fc=i Z7r 



is conserved after this transformation. 

Because V'm.m ^ s a ^ so an eigenfunction of the integrals of motion A and A* with their 
eigenvalues \ m and AJ„ = A~: 

^m,m = A m^ m ,m> -4* ^m.m = A m ^m,m > A = Pl2-PnlPl-Pn , (77) 

one can verify that, for the unitarity of the duality transformation, the constant c ~ should 
be chosen as 

C m,m = l^m| 2™ , (78) 

for an appropriate phase of ip m ^. Here the factor 2™ appears due to the relation d 2 „ = Add*. 
This value of c m ~ is compatible also with the requirement that two subsequent duality 
transformations are equivalent to the cyclic permutation % — > % + 1 of gluon indices. 

Thus, the duality symmetry takes the form of the following integral equation for ip m ~: 



f n ~ 1 d 2 p' k _ lk n e^fc.fc+iP'fc -y> > 

i> m ,miPl2,---,p2) = \\m\ 2" / J] 0qT ' II -j 2-1>k,m(Pw> ->Pnl)- ( 7 9) 

*=i ^ *=i PU+i 
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Note that the validity of this equation in the Pomeron case n = 2 can be verified from 
the relations 



/m,m(Pl) P^'i A)) 



P12 \ ( P*12 \ 
P10P20/ VPI0P20/ 



I A m ,| 2' 



/ 



2tt 



1/012 1 



(P / I 2 ) m (p'i2) m = e^ m )pi2 1 - m A^ 1 " m , M = - m)| , 



i8(m,m) 



- 2 m+m-i f f m*(m*-l) \ 
1 ; r(2 - m) V - 1) 7 



1/2 



(80) 



Let us use for f m ~ the conformally covariant anzatz 



/m,m(Pl,-,At;Po) = LQ 



\ m/n / n „* \ «V n 
\ / TT Pi,i+1 | 



n 



. I. PiO 



/m,m(^l, -,^) , ( 81 ) 



\i=l ^iO 

where the anharmonic ratios x r (r = 1, 2...n) of the gluon coordinates are chosen as follows 

Pr-l,r Pr+1,0 



x r = - ; !I-'V ( D": E(-!) r 11 *k = 0- 



Pr-lfiPr+l,r r= \ r =l fc=r+l 

The function f m ~(xi, ...,xt) is invariant under certain modular transformations as a conse- 
quence of the Bose symmetry. 

For the physically interesting case m, m — > 1/2 we can calculate tp m ^, taking into account 
the logarithmic divergence of the integral at po - ► 00: 



m,m->l/2 1 — m — m 

Let us change the variables z r to new ones as follows: 

n— 1 n 

\n— fe— 1 



Pr-l,r 
Pr+1, 



(83) 



J/k = Pk+l,k/Pn,n-l = (-1)™ fe 1 LI *r > J/n-1 = 1, S J/* = 0. (84) 

r=fe+l fe=l 

The duality equation for the wave function f(yt, ...,y n -2) can be written in the form: 



. n-2 j 2 „,/ 



il/n 



J fc=i ^ k=i ml 

The integral kernel K(lf; y ) is given below: 

K(f,7) = rJ 



(85) 







n-1 






3 


2tt 


Pn,n—1 


\Pn,n—l \ 



— exp ^-i £ (pi, y fe -lPn,n-l)) 



(86) 



and is calculated analytically: 



Vk-iPk 

k=i 



In l _1 — _9 T 

\l J n,n— 1| — ^ 



n-1 fc-1 
fe=l r=n 



(87) 



To simplify the duality equation, we consider below the compound state of three reggeized 
gluons. 
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3 Duality equation for the Odderon wave function 

In the case of the Odderon the conformal invariance fixes the solution of the Schrodinger 
equation [5] 

(\ m/3 / * * * \ m/3 

Pl2P23P3l\ P12P23P31\ , /->n /ooN 

n 2 n 2 n 2 n *2 n *2 n *2 /m,m( X J l 8S J 
Pl0P20P30/ \F10F20F30/ 

up to an arbitrary function / mm ("a^) of one complex variable x being the anharmonic ratio 
of four coordinates 

x= P -^. (89) 
P10 P32 

Note that, owing to the Bose symmetry of the Odderon wave function, / mm (a^) has the 
following modular properties: 

fmd^) = (-1) ( " )/3 fmd^/ l^l 2 ) = (-1) ( " )/3 - *) (9°) 

and satisfies the normalization condition 

f ~ 2 = f tl f ~<-i) 2 (91) 

Jm,m J 1 , \| 4 / 3 Jm < m ^ ' ' \ U± J 

compatible with the modular symmetry. 

After changing the integration variable from p to x in accordance with the relations 

P3lPl2 , P3lPl2P32 7 / no x 

Po = h/Oi, dp = -, da; (92) 

Pi2-a:p32 \P\2-xpm) 

the Odderon wave function i/) m m(f)ij) at if = can be written as 

Til — 1 / ^ \ TO— 1 



/ \ Til — 1 / „. \ 1 



where 



2n\x-z\* \x(l-x)J \x*(l-x*) 
In fact this function is proportional to f\- m \-zJC^)'- 

Ut) - (x(l - x)) 2 ^ 1 ^ 3 (*•(! - x*)) 2 ^/ 3 A_ m)1 _ m (y) , (95) 

which is a realization of the discussed linear dependence between two representations (to, to) 
and (1 — to, 1 — to). 

The corresponding reality property for the Mobius group representations can be presented 
in the form of the integral relation 

Xto,to(^) = / ^ (* - ^"^ (*• - ^) 2 ™" 2 Xi-to,i-to(^) (96) 
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for an appropriate choice of phases of the functions X m ,m an d Xi- m ,i-m- The function 
Xm,m(~^) satisfies the modular relations 

Xm,m( Z ) = ( — 1) Z Z Xm,m( z I \ z \ ) = ( — -Q Xm,m( ^ ~ Z ) • (97) 

The duality property of the wave function ^ m ^{Pij) can be written in a form of the 
integral equation: 

-fTT^-IX I 9 3 / rf Vl2 ^23 eX P(^(/^ Pg + /4 P3l)) /T/N fqoN 

- |A m | 2 y -^-^ ^r^rp > ( 98 ) 

where |A m | 2 is the corresponding eigenvalue of the differential operator 

|-4| 2 = \pl2p23p3lPlP2P3\ 2 • (99) 

The function ^~ >m (p]J) is transformed as ^i-™,i- m (pj): 

^(Pj) ~ ^l-m,l-m(^) • ( 10 °) 

In terms of Xm,m(^) ^ ne above duality equation looks as follows 

w ™-/ ™rlPT?£ «*- 7 *- {7) - (101) 

J 27r(z(l — z)J (z*(1 — z*J) 



|A m | 



where 

R( ~z )= I 8d / P f (^ 2l $ (p32p3 2 ) m {PI2P32T exp(i(pi2p' 12 + Pnp'13)) • ( 102 ) 

J |P32l I P32 I 

This integral is calculated to be 



with 



R(l?,7) = C m ~Z l - m (Z*) 1 -™ (103) 

^ 2e a (m ,£) ^ (_gg-rn r(1/2 ) r(l/2) 

m ' m |m(l-m)| sin(Trm) T(2 - m)T(2 -m) ' 1 J 

where the phase 5(m,m) was defined above and 

Z = zz'* -z+l = z(z'* - 1 + i) . (105) 

Further, we introduce the new integration variable z 1 — > z'*, which effectively leads to the 

substitution Xm.m( z ) ~^ X m m( z )■ By performing the modular transformation 1 — ^ — > z 
and taking into account the modular relation 

Xm^) = (1 " ^^^(l - O^UO* " *)/ I 1 " ^l 2 ) > ( 106 ) 
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one can rewrite the above duality equation for X m ,m(~^) as 



X 



m,m \ 



IX C ~ 



d 2 z' (z'(\-z'\z(\-z)\ m - X fz'(l-z')z(l-z)Y m ~ 1 



2n 



z' — z 



z' — z 



X*d~2)- (107) 



This equation for x m m (^) corresponds to the symmetry of the Odderon wave function under 
the involution p k \skimPim- 

It can be written in the pseudo-differential form: 



z(i - z) (idy-™ z *(i - z *) wy- m Vl _ mtl _~o?) 

where 



A 



m,m 



(vW-mC*))' , (108) 



^i-Sl *>) = 2 (1 - m - m)/2 (*(1 " ^)) 1 " m " ^*)) i " m ^s(t) • (109) 

Note, that for a self-consistency of the above equation its right-hand side should be ortogonal 
to the zero modes of the operator (id) 2 ~ m (id*) 2 ~ m . Due to the Bose symmetry of the wave 
function it is enough to impose on it only one integral constraint: 



rn 



I \z(l-z)\ 2 ( V?1 —' 1 - 



:(l-z)\ 2 

The normalization condition for the wave function 

2 r d 2 x 



0. 



(110) 



\x(l — x) 



is compatible with the duality symmetry. 

The holomorphic and anti-holomorphic factors of / mm (pi, Pi, ~f>z\ Po) are eigenfunctions 
of the integrals of motion A and A*: 

Af m = \ m f m , A*f~=\~f~, \~ = (\ m )*. (in) 

In the limit m — * 1/2, m — ► 1/2, corresponding to the ground state, the integral over x 
for the wave function i/) m ~ at q = is calculated, since the main contribution appears from 
the singularity at x = z: 



^m,m(P?j) - 



1 /3 

Pl2P23P3l| f{Z,Z) , Z = p 12 /p 3 2, 



(112) 



(m + m — 1) 

where f(x,x*) = f\/2,i/i{x,x*). This means that one can obtain for f(x,x*) the following 
equation in the x-representation: 

| 1/3 f(x,x' 



\x(l — x) 



4 |A| r d 2 y \y - x\ 



* J \v(l 



f(y,y* 



(113) 



W-y)f /3 

for a definite choice of the phase of the ground-state function /. 

The common factor in front of the integral is in agreement with the condition of conser- 
vation of the norm of the wave function after this canonical transformation (only its sign 
can be changed). The above integral relation is reduced to the following pseudo-differential 
equation 
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2 

x(l-x)p 3/2 ip(x,x*) = \\\ ip(x,x*), (114) 
where p = id/d(x) , p* = id/d(x*), if we introduce the new function 

ip(x,x*) = \x(l -x)\ 1/3 f{x,x*) (115) 

and use the relation 

\P\ m = -\p\ / 7T- / — ^-y exp(— 2 p X ) — — — o (x ) . (116) 

J-ir/2 271 J-oo \p\ 4 

For self-consistency of the pseudo-differential equation one should impose on the wave func- 
tion the following constraint: 

<Hr cp(x,x*) = 0. (117) 



\x(l — X 

Let us derive the duality equation for (p mtm (J?) f° r general values of m and m by using 
different arguments. We start with the conformally covariant anzatz for the holomorphic 
factor / m (pi,p2,P3;po): 

(\ m/3 
P12P23P31 \ , / \ P12P30 / 110 \ 
~2 2 T fm{x) , X = . (118) 

PlO P20 P30/ PlO P32 

In the x representation the integral of motion A = p 12 P23 P31 P1P2P3 is an ordinary differ- 
ential operator. When acting on f m (x), it can be presented in the following form 



- (-(x - 2) + x(l - x)d) (-(1 + x) + x(l - x)d) - (-(1 - 2x) + x(l - x)d) 

x V 3 / X x V 3 / x V 3 / 

= X-™/ 3 A m X m / 3 , X = x (1 - x) , A m = ai . m a m , (119) 

where 

a m = x (1 -x)p m+1 , a^ m = x (1 - x) p 2 ~ m , p = id. (120) 

Therefore the differential equation A f m = \ rn f m for eigenfunctions f m and eigenvalues 
\ m is equivalent to the system of two pseudo-differential equations 

dm <Pm = Kn fl-m i a l-m <Pl-m = h-m , (121) 

where we introduced the functions ip m and <pi_ m in accordance with the definitions 

Vm = X m ' 3 f mi vi-m = X^ m ^ A_ m (122) 
and the eigenvalues l m and l\- m are related with the eigenvalue \ m of A m as follows 

ll—m,lm = \n- (123) 

The function ipi_ m (x) has the conformal weight equal to 1 — m. It is important that 
there is another relation between ip m and y?i_ m 

(fl-m = R m (x,p) (fm , (124) 
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where 

R m (x, P ) = x- m+l P l ~ 2m X~ m , i2x_ m (x,p) = (Rm(^P)r 1 ■ (125) 

This follows from the fact that, for the Mobius group, the complex-conjugated representa- 
tions O m ~(po) and O 1 _ ml _~(po) are linearly dependent. To obtain the above expression 
for R m (x,p) one should use for the Odderon wave function f m ^(p~i, pi, pi; Po) the coformal 
anzatz and the property of holomorphic factorization. The transformation 

<p'(l>) = U<p(l>) , U = R m (x,p) R~(x*,p*) (126) 
is unitary for our choice of the norm: 



|<p'U 2 = IMI 2 ■ (127) 



Because 

d 2 x 



— (<p(^))* A m (x) A~(x*) <p(&) = \\A m <p\\ , 

x(l — x)\ 

the eigenvalue of the operator A~(x*) is complex conjugated to A m : 

A m (x) <p(-&) = X m <p(-&) , A~(x*) <p(-&) = X* m <p(-&) . (128) 

Note that due to its Mobius covariance (p(l?) is the sum of products of the eigenf unctions 
having opposite signs of their eigenvalues A: 

= E C * (<a(z) vi,A X *) + <Pm,-x( X ) ^.-A'OO ) • 

Since under the simultaneous transformations x <-> re* and m <-> m this function should be 
symmetric for fixed A, we conclude, that eigenvalues satisfy one of two relations 

A m = ±(A m )* (129) 

and therefore A m can be purely real or imaginary. It turns out, that A m is imaginary as a 
consequence of the modular invariance [9] . 

One can veryfy the validity of the following relation 

Rm(x,p) Ql—m dm &m 0,1— m Rm(x,p), (130) 

if the following identity is used: 

X~ m ai _ m a m X m = x 2 (l- x) 2 p 3 + 2(1 + m) x(l -x){ip{\- 2x) p - p) + 

m(l + m) (x(l — x)p — (l-2x) 2 p + im(l-2x)Y (131) 
In particular this relation means that two eigenvalues of A m coincide 

A m = Ai_ m (132) 

for the eigenfunctions ip m and <pi_ m linearly related by R m (x,p). Let us introduce the 
operators 

S m (x, p) = Ri- m (x, p) a m , S[_ m (x, p) = ai_ rn R m (x, p) , (133) 
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leaving the function ip m in the same space. Due to the above formulas these operators 
commute with one another: 

A m = Sl_ m {x,p) S m {x,p) = S m {x,p) Sl_ m (x,p) . (134) 

Therefore we obtain two duality equations for the wave function whose conformal weight 
is equal to m: 

S m {x,p) if m = (fi m , S[_ m (x,p) if m = L { ™ ] ip m . (135) 

As a consequence of the relation 

L (m) L (m) = A (136) 

between the eigenvalues of S m and S{_ m , the second equation follows from the first one. 

In the particular case m = 1/2, the equation for the holomorphic factors <p(x) looks 
especially simple: 

x(l - x)p^ V±v - x {x) = ±V\cp ±v - x (x) (137) 

and can be reduced in the p-representation to the Schrodinger equation with the potential 
V{p) = p^l 2 . For m = m = 1/2 the total odderon wave function ip(x, x*), symmetric under 
the above canonical transformation, is a solution of the equation: 

\x(l-x)\ 2 \p\ 3 <p(x,x*) = |A| <p(x,x*), (138) 

where in accordance with its hermicity properties the eigenvalue of the operator y^A* for the 
anti- holomorphic factor (p(x*) is taken to be equal to A*. 

4 Single-valuedness condition 

There are three independent solutions ip\ m \x, A) of the third-order ordinary differential equa- 
tion 

Gi-m flm f — — x) (x(l — x)d 2 + (2 — m) ((1 — 2x)d — 1 + to)) d (f = A ip (139) 

for each eigenvalue A. In the region x — > they can be chosen as follows 

oo 

4 m) (^ A) = £ 4 m) (A) x k , 4 m) (A) = 1 . (140) 
k=i 

oo 

^ m) (^A) = ^4 m) (A)x fc +^(x,A)lnx, = , (141) 

k=0 

oo 

A) = £ cKJA) ^ , 4 m) (A) = 1 • (142) 

k=0 

The appearance of In a; in <^( m )(x, A) is related with the degeneracy of the differential equation 
in the small-x region. There is an ambiguity in the definition of ip( m \x, A) because one can 
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add to it the function (p( m \x, A) with an arbitrary coefficient. We have chosen a™ = to 
remove this uncertainty. 

Taking into account that for p 1 — > p 2 the operator product expansion is applicable, the 
functions tp\ m ^ (x, A) can be considered as contributions of the holomorphic composite opera- 
tors 0\ M \pi) with the conformal weights M = 0, m and 1 for i — s, f and r correspondingly. 
In this interpretation the above degeneracy is related with the existence of the conserved 
vector current (for m — 1/2 there is also a conserved fermion current). 

Due to the above differential equation the coefficients a, c and d satisfy the following 
recurrence relations 

iXa^ = (4:1 + A) 

+ (4-1 + 4-1^) (k - l)(k -m)(k-l-m), 

iXc^l = (k + m)(k + m + l)(k + l)c£l +1 -(k + m)k(2k + m)c { ™l 



+ (k + m- l)k(k- l)c (m) 



k+m 

k+m—1 ' 



i\d { ™ } = k(k + l)(k + l-m)d ! il\-k(k-m)(2k-m)dt ) (143) 
+(k-l)(k-m)(k-l-m)d { ™\. (144) 

In particular, from the equation for 4"^ & t k = 0, since 4"^ = 1, we obtain 

4 m ) = L ( m _ i) . (145) 
A 

The introduced functions have simple analytic properties in the vicinity of the point 
x = 0. In particular, A) is regular here and is transformed under the modular 

transformation 

x -> x' = -x/{l - x) (146) 

as 

ipl m \x', A) = -(1 - x) m ipi m \x, -A) . (147) 

The functions ^^(x, A) and ip^ n \x, A) have singularities at x — 0, which leads to differ- 
ent results for their analytic continuations to negative values of x: 

^("V, A) = -(1 - x) m -A) ± m 4 m \x, - A)) , 

i P i f{x',\) =exp(±i7rm)(l-x) m </4 m) 0z;,-A). (148) 
Therefore, from the Bose symmetry of the Odderon wave function 



rn 



/m,m(Pl»P2»P3;Po) = ( I ( f ^ ) ^^l, ( 149 ) 

\ P20 P30 / \ P20 P30 / 
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combined with the single- valuedness condition near ~af = 0, we obtain for the total wave 
function the following representation: 

U 1 - **) = *T\ X , A) <pf\x\ A*) + d A) A*) + ^(s, A) A*)) 

+ c 2 ^(x, A) 4™\x*, A*) + (A -> -A) . (150) 

The complex coefficients ci, c 2 and the eigenvalues A are fixed from the conditions of the 
single-valuedness of / m ^(pi, P2, P3; Po) at pjj = pi (i — 1,2) and the Bose symmetry [9]. It 
is sufficient to require its invariance under the transformation p% <-> p| corresponding to the 
symmetry of </? m ~(a;, x*) 

( Pm,m( x > x *) = <Pm,m( 1 - x > 1 -X*). (151) 

For this purpose, one should analytically continue the functions ip\ m \x) in the region near 
the point x = 1 and calculate from the differential equation the monodromy matrix 
defined by the relations 

^ m) (*,A) = EC^ri m) (l-^-A), (152) 
k 

#)(**, A*) = EcfVPa-^-A*). (153) 

Owing to the single-valuedness condition and the Bose symmetry of f m ~ we obtain a set 

of linear equations for parameters C\ and c 2 with coefficients expressed in terms of and 

Cl™\ The spectrum of A is fixed by the condition of selfconsistency of these equations [9]. 

To derive the relations among parameters c±, c 2 and A following from the duality symme- 
try, it is convenient to introduce the operators 

S m fn = S m (x, P ) S~(x*, P *) , ~ = S[_ m {x, P ) S[_~{x\ P *) , (154) 

where S m (x,p) and S{_ m (x, P ) were defined in the previous section. These operators are 
hermitially conjugated each to another and have the property: 

\A m (x)\ 2 = S+ ~ S m ~ = S m ~ S+ ~ . (155) 

In particular it means, that they have common eigenfunctions 

S m ,m <Pm,m( X i X *) = ^ <Pm,m( x , x *) , S+~ (f m ^(x, X*) = C X e~ l6 <P m ^(x, X*) , (156) 

where 

ao _ ( _^m-m T(m) r(l + m)r(l + m) 

~ [ ] r(l-m)r(2-m)r(2-m)- 1 ' 

The eigenvalues are obtained from the small-rr asymptotics of these equations. 
Because the operators S m ~ and 5+ ~ are hermitially conjugated, we have 

|ci| = |A|. (158) 
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In the particular case m = m = 1/2, where Si/2,1/2 = | -S'i /2 1 2 , the coefficient C\ is positive: 

ci = |A| . (159) 

Another relation can be derived if we shall take into account, that the complex conjugated 
representations ip m ~ and <p 1 _ m x _~ of the Mobius group are related by the unitary operator 
U = R m (x, p)R~(x*, p*), defined in the previous section: 

e ll (v m ^r = U Vm ~, (160) 

where e n is an eigenvalue of this operator. By calculating the right hand side of this equation 
at x — > 0, we obtain 

e- = (-1)— r( 2 -rn)r(2-m) = 

v ; r(i + m)r(i + m) c\ ' ci v ; v ; 

One can verify from the numerical results of ref. [9] that both relations for c\ and c 2 are 
fulfilled. For example, we have for the ground-state eigenfunction with m — m = 1/2: 

i A = 0.205257506 , a = 0.205257506 (162) 

and for one of the excited states with m — m = 1/2 + % 3/10: 

% A = 0.247227544 , Cl = 0.247186043 - i 0.004529717 , c 2 = -1.156524786 - i 0.415163678, 

|ci| = |A| , Im— = 27m (1/2 + iS/10)' 1 . (163) 

Cl 

After the Fourier transformation of y^m.ml"^) to the momentum space the regular terms 
near the points "of = and it = 1 do not give any contribution to its asymptotic behaviour 
at — > 00. The requirement of the holomorphic factorization and single- valuedness of the 
wave function in the momentum space leads to the quantization of A. 

We can obtain from the duality equation and reality condition also the representations 
for the coefficients Ci^ in terms of integrals from {p m ~(x, x*) over the fundamental region 
of the modular group, where the expansion in x is convergent. These relations allow one to 
calculate the coefficients Ci j2 without using the single-valuedness condition. 

5 Hamiltonian and integrals of motion 

The holomorphic Hamiltonian h for the compound state of n Reggeons for iV c — > 00 com- 
mutes with the transfer matrix T{u) owing to the following relation for h^k+i'- 

[h ktk+1 ,T(u)] = -itr (L 1 (u)...L k _ 1 (u) (L k (u) - L k+1 (u)) L k+2 (u)...L n {u)) . (164) 

It can be considered as a linear equation for h k ,k+\- The formal solution of this equation can 
be written as 

h ktk+1 =lim(ij*dt' exp{tT{u)t')[h ktk+1 ,T{u)]exp{-tT{u)t') + h k , k+1 {t)^ , (165) 
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where h k>k+ i(t) is the time-dependent operator 

h k ,k+i{t) = exp(iT(u)t) h k:k+1 exp(-iT(u)t) . (166) 
Since the integral term is cancelled in the sum of h ktk+ i, we can substitute 



fc+i 



h 



k,k+l 



(t). 



(167) 



At t — > oo as a result of rapid oscillations of off-diagonal matrix elements, each pair Hamil- 
tonian is diagonalized in the representation, where the transfer matrix is diagonal, and 
therefore it is a function of the integrals of motion q^: 



h kyk+1 (oo) = fk,k+i(q2,q3, -q n ) 



(168) 



Its dependence from the spectral parameter u disappears in this limit and the total Hamil- 
tonian is 



h = h(q 2 , q 3 , ...q n ) = /k,*+i(&» qs, ...q n ). 



k=i 



All operators 0(t) satisfy the Heisenberg equations 

-ij-0{t) = [T{u),0{t)] 



(169) 



(170) 



with certain initial conditions. In the case of the pair Hamiltonian the initial conditions are 

/ifc,fc+i(0) = ip(m kjk+1 ) + ^(1 - m fc>fc+ i) - 2^(1) , (171) 
where the quantities rfik,k+i are related to the pair Casimir operators as 

m Kk+1 (m k , k+1 - 1) = Ml k+l = -p 2 kjk+1 d k d k+1 . (172) 
In the case of the Odderon, h does not depend on time if h k>k+ i{t) is determined as 

h k , k+1 (t) = e UA h k , k+1 e- UA , (173) 

and h k)k+ i(oo) is a function of the total conformal momentum M = m(m — 1) and of the 
integral of motion q 3 = A, which can be written as follows: 



A = - 
2 



M 2 12 , M 2 13 



Mi 3 , Ml 2 



M?3 , M 2 2 3 



(174) 



Using these formulas and the following relations among the Mobius group generators ~E% r 



Ml - Ml 



hi kl 



2 (A?* - Ii k , 
Ml, M t -M k }} =A(M % -M k ) , 



(175) 
(176) 



we can verify the commutation relations 
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% [h 12 , A] 



M 2 3 - M| 3 , i [h 13 , A] = M| 3 - M 2 2 , % [h 23 , A] 



M 2 2 



(177) 



from which it is obvious that A commutes with h. 

In a general case of n reggeized gluons, one can use the Clebsch-Gordan approach, based 
on the construction of common eigenfunctions of the total momentum M and a set jMfcj 
of the commuting sub-momenta, to find all operators M kk+1 in the corresponding repre- 
sentation. However to calculate h we should perform an unitary transformation to the 
representation, where T(u) is diagonal, because in this case for t — > oo the off-diagonal 
matrix elements of M kk+1 disappear and their diagonal elements depend only on q r : 



fk,k+l(<l2, 93, — Qn) — (?2, ■■■Qn | ^fc,fc+l | 92, ■■■Qn) , 



Qk | 92, ■■■Qn) = Qk | 92, ...g n 



(178) 



Let us consider, for example, the interaction between particles 1 and 2. The transfer 
matrix, which should be diagonalized, can be written as follows 



T{u)= (u 2 - - L )d 3 ... n (u) + (iu L - - 



) ^3...n(«) , 



(179) 



where the differential operators d$_ n (u) and d 3 ..„(«) are independent of pi and p|. They 
are related to the monodromy matrix t 3 „ n {u) for particles 3,4, ...,n as follows 



d 3 ... n {u) = trt 3 . 



u 



, u 3...n 



£r( a%... n (u)), t 3 ...„(u) = L 3 (u)...L n (u) 



(180) 



and the matrix t 3 „ n (u) satisfies the Yang-Baxter equations with a hidden Lorentz symmetry. 

The operators L and are constructed in terms of the Mobius group generators of 
particles 1 and 2: 



? = tf 1 + jZ 2 , 7$ = ^-^, M z k =p k d k ,M k + = -p 2 k d k ,M k =d k . (181) 
They have the commutation relations, corresponding to the Lorentz algebra: 



L 2 , L ± 



2iV 2 



L+L" 



N z , N ± 



2L d 



±L 



± 



7V+ AT 



±iV 



± 



2L 2 



(182) 



Let us introduce the Polyakov basis for the wave function of the composite state of two 
gluons with the conformal weight M: 



Pv,M) = 



Pl2 



M 



(183) 



^PlO'P20' ) 

Here po' enumerates the components of the infinite-dimensional irreducible representation of 
the conformal group. 
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One can verify that the representation of the generators L and in this basis is given 

by 

M z | p Q/ ,M) = ( Pl d 1 + p 2 d 2 ) \p ,,M) = -(p ,d , + M) \pv,M), 
M + | pv, M) = -(pjd, + p\d 2 ) | pv, M) = (p 2 0l d , + 2Mpff) | p<y, M), 

M~ | p >, M) = {8 1 + d 2 ) | p<y, M) = -d ' | po> , M) (184) 

and 

(N z -p Q ,N-) I Ay ,M) = -^ T ^| Ay ,M-l), 
(iV + + 2p o; iV z -p2 ;A r-) | p ,,M) = -2M|p o; ,M-l), 

2M - 1 1 

iV" \pv,M)=\pv,M+l)+ — n . 2 4 |po',M-l). (185) 



M(M-l) Iru ' ' Iru ' ' (M-l) 2 

Note that there is a simple relation among the generators, provided that they act on the 
state | po>, M): 

N z - pv N~ , N + + 2p ' N z - p 2 ,N~] = M + + 2p , M z - p^AT = . (186) 

The eigenfunction of the transfer matrix T(u) can be written as a superposition of the 
states | po',M) with various values of po> and M: 

fm(pl,p 2 ,...,p n ;Po)=J2 / d P0' \P0',M) f m , M (pO',p3, ■■; Pn, Po) , (187) 

M J 

where m is the conformal weight of the composite state. The function f m> M(po'---Pn', Po) i n 
accordance with the Mobius symmetry, has the form 

m + M 

M 1 n ( PrO \ ™~ 2 

/m,M(P0',P3,--,Pn;P0) = (P 'o) _m+M_1 II ~ ^1, X 2 , X n _ 3 ) , (188) 



r=3 



,PrO' 



where x r are independent anharmonic ratios constructed from the coordinates Po',P3, ■■■,Pn- 
Because of its Mobius invariance, the transfer matrix T{u) after acting on / m (pi, p„; po) 
gives again a superposition of the states | p ', M), but with the coefficients which are linear 
combinations of f m ,M and / m ,M±i- Therefore for its eigen function the coefficients satisfy 
some recurrence relations, and the problem of the diagonalization of the transfer matrix T{u) 
is reduced to the solution of these recurrence relations. For n > 3 in the sub-channel p\ t2 
the recurrence relations depend on matrix elements of the operators a^„_ n {u) and ds,„ n (u) 
between the wave functions f m ,M(po', P3, Pn', Po) which should be chosen in such a way, 
to provide the property of / m (pi, p 2 , ■■■,p n ] Po) to be a representation of the cyclic group of 
transformations i — > i + 1. 

In the appendix we consider these relations in the first non-trivial case n — 3. In the 
next section the relation between the Odderon Hamiltonian and its integral of motion A is 
discussed from another point of view. 
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6 Odderon Hamiltonian in the normal order 

Let us write down the pair Hamiltonian as follows [4] 

h 12 = \og(p 2 l2 d,) + \og(p 2 l2 d 2 ) - 2 log( Pl2 ) - 2^(1) . (189) 

This representation allows us to present the total Hamiltonian for n reggeized gluons in the 
form invariant under the Mobius transformations 

h=±Lg( Pfc+2 '° Plk+1 d k ) + log ( Pk - 2 >° d k ) - 2 m) , (190) 

k=l \ \Pk+l,0 Pk+l,k+2 ) \Pk~l,0 Pk-l,k-2 ) ) 

where po is the coordinate of the composite state. 

We consider below in more detail the Odderon. Using for its wave function the conformal 
anzatz 

— — j (p m (x), x = 12 30 , (191) 
P20P30 / P10P32 

one can obtain the following Hamiltonian for the function (p m (x) in the space of the anhar- 
monic ratio x [4] 

h = 67 + log (x 2 d) + log ((1 - xfd) + log - x ) d + + 

log - x ) d + m ))) + lo § ( - ~^~ ( xd ~ + log {l^ xd ~ m • ( 192 ) 

It is convenient to introduce the logarithmic derivative P = xd as a new momentum. 
Using the relations [4]: 

log{x 2 d) = \og(x) + - P) , log(<9) = - log(x) + i/>(-P) , 
log(x 2 d) = log(9) + 21og(a;)-^, 

((1 -x)d + m) = (1 -x) 1+m d{\ -x)- m , xd - m = x 1+m d x~ m , (193) 
one can transform this Hamiltonian to the normal form: 

h 00 

_ log^) + ^(1 _ P ) + + _ P ) _ 3^(1) + ^x k f k (P) , (194) 

1 k=i 



where 



r /r,N 2 1/ 1 1 \ * Ct(fc) , N 

a ^ = -* + 2(pTI^ + pTi) + SfT7' (195) 
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Here 

_ (-1)*-* r(m + t) ((t - fc) (m + t) + m fc/2) 

H J fcr(m-ib + t + i)r(t + i)r(A;-t + i) ' 1 J 

On the other hand the differential operators a m and ai_ m can be written in terms of the 
quantities P and x as follows 

a m = i-^x-^l-x) ^-^^ , (197) 

ai-m = r 2+ -x- 1+m (l-x) r( ~^~^ +2) . (198) 

Using the above representation for h and the following expression for the integral of 
motion: 

Cl — 

P = % ai_ m a m = - - ((1- x)P -1- x + xm) P (P - m) , (199) 

one can verify their commutativity 

[h, B] = 0. (200) 

Therefore h is a function of B. 

In particular for large B this function should have the form: 

\= log(S) + 3 7 + £ . (201) 

The first two terms of this asymptotic expansion were calculated in ref. [4]. The series 
is constructed in inverse powers of B 2 , because h should be invariant under all modular 
transformations, including the inversion x — > 1/x under which B changes its sign. The same 
functional relation should be valid for the eigenvalues e/2 and /i — % A of these operators: 

£ 00 r 

- = \og(n) + 3 7 + £ ;i • (202) 

Z r=l P 

For large \x it is convenient to consider the corresponding eigenvalue equations in the P 
representation, where x is the shift operator 

X = 6XP( ~^ } ' (203) 
after extracting from eigenfunctions of B and h the common factor 

ip m {P) = T(-P) r(l - P) r(m - P) exp(mP) $ m (P). (204) 

The function 3> m (P) can be expanded in series over \j\x 

oo 

$ m (P) = E //-"^(P) , C(P) = 1 , (205) 

?1=0 
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where the coefficients $^(P) turn out to be the polynomials of order 4n satisfying the 
recurrence relation: 

®!{P) = £(* - 1) (* - 1 - m) ((* - m ) $ ™ ^ - 1) + (* - 2) - 1 - m)) 

fc=i 

- - - 1) (k - 1 - m) ((k - m) ^ m -\k - 1) + (k - 2) ^(fc - 1 - m)) , (206) 
z fe=i 

valid due to the duality equation written below for a definite choice of the phase of $ m (P) 

aT m (1 - x P (P - m) (P - m + 1)) $ m (P) = //" $!_ m (P) (207) 

with the use of the substitution x \i — > x. 

Note that the summation constants $„(0) in the above recurrence relation have the 
anti-symmetry property 

C(0) = ?_ m (0) , (208) 
which guarantees the fulfilment of the relation 

*S,M = *i-m(°) (209) 
being a consequence of the duality relation 

ct^(P) = ^_ m (P - m ) + (P-l)(P-m)(P-m- l)^\P - 1) . (210) 

The symmetric part of $^(0) under the substitution m <-> 1 — m would simply modify the 
normalization constant for $ m (P). 

The energy can be expressed in terms of $ m (P) as follows: 

£ -= log( / i) + 3 7 + ^plog$ m (P) 

oo k 

+ (^m(P)r 1 E UP - k) $ m (P - k) H(P - r){P - r + 1)(P - r - m + 1) (211) 

k=l r=l 

and it should not depend on P due to the commutativity of h and B. 

By solving the recurrence relations for $„(P) and putting the result in the above ex- 
pression, we obtain the following asymptotic expansion for e/2: 

£ = , og( ,) + 37 + (A + J£ (m _ 1/2) . _ _L (m _ 1/2) «) i + 



/ 4185 2151 . , /0 , 2 \ 1 / 965925 \ 1 

("2050048 - 49280 (m " 1/2) + -) J? + 137044224 + "J J? + " ' (212) 

This expansion can be used with a certain accuracy even for the smallest eigenvalue \x = 
0.20526, corresponding to the ground-state energy e = 0.49434 [9]. For the first excited 
state with the same conformal weight m = 1/2, where e = 5.16930 and ji = 2.34392 [9], 
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the energy can be calculated from the above asymptotic series with a good precision. The 
analytic approach, developed in this section, should be compared with the method based on 
the Baxter equation [10]. 

In the conclusion, we note that the remarkable properties of the Reggeon dynamics are 
presumably related with supersymmetry. In the continuum limit n — > oo the above du- 
ality transformation coincides with the supersymmetric translation, which is presumably 
connected with the observation [11], that in this limit the underlying model is a twisted 
N = 2 supersymmetric topological field theory. Additional arguments supporting the su- 
persymmetric nature of the integrability of the reggeon dynamics were given in ref. [12]. 
Namely, the eigenvalues of the integral kernels in the evolution equations for quasi-partonic 
operators in the N = 4 supersymmetric Yang-Mills theory are proportional to — 1), 
which means that these evolution equations in the multicolour limit are equivalent to the 
Schrodinger equation for the integrable Heisenberg spin model similar to the one found in 
the Regge limit [7]. Note that at large N c the N = 4 Yang-Mills theory is guessed to be 
related with the low-energy asymptotics of a superstring model [13]. 
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Here we consider consequences of the conformal weight representation for the Odderon. In 
this case the total Hamiltonian is 



Appendix 



h(m, A) = hi2 + h 2 3 + h 
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The eigenvalue of h is expressed in terms of its matrix elements: 



3 



h(m, A) = ^2(m, A | h k ,k+i \m,X), 



where | m, A) is a normalized eigenfunction of two commuting operators 




m, A) 



m(m — 1) | m, A) , A \ m, A) = A | m, A) . 
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Let us consider for deflniteness the interaction in the channel 12, where M is the pair 
conformal weight. If one will construct the matrix V^(m) performing the unitary transfor- 
mation between the M- and A-representations 

\m,\) = Y,V&(m)\m,M), M 2 2 | m, M) = M(M - 1) | m, M) , £^l^ = ^„ 

M M 

then the diagonal matrix elements of the pair Hamiltonian h i2 can be calculated as 
(m, A | h 12 | m, A) = ^ h(M) (m)V^(m) , h(M) = ip(M) + ^(1 - M) + 2 7 . 

M 

We shall derive below the recurrence relations for V^. To begin with, we note that, 
according to the commutation relations for the Lorentz algebra generators, there are non- 
vanishing matrix elements of the boost operator only between the states | M) and | Mil). 
This is valid also for the matrix elements of the operators A and M 2 3 — Mf 3 according to 
the relations 

(M'\A\M) = -(M 1 -M)(M' + M- 1) (M' | M 2 3 - M 2 3 | M) , 

(M ± 1 | M 2 3 - M 2 3 | M) = 2Ml (M ± 1 | jtf | M). 
Thus, for the common eigenfunctions | to, M) of two commuting Casimir operators 

(^Mk^j | m,M) =m(m- 1) | m,M) 

and 

M 2 2 \m,M) = M(M - 1) | m, M) 

we have 

A | to, M) = t (MC+ M | to, M + 1) - (M - 1)C- M | m, M - 1)) , 
where the coefficients C^ M are defined by the relations 

(M 2 3 - M 2 2 3 ) | m, M) = C+ M | m, M + 1) + C mjM | m,M — 1) . 

One can obtain from the above equations the following recurrence relation for the unitary 
matrix V^(m): 




To calculate the matrix elements C^ M of the operator M 2 3 — M| 3 , we use the above 
representation of the generators L and N in the Polyakov basis and obtain 

(M 2 3 - M 2 2 3 ) | p ', M) = (27V 2 M| + N + M~ + N~M^ \ p >,M) = 

27 



2M 



— j p3o>d > -Ij | p<y, M - 1) d 3 - p 2 w N I p ', M) d 3 



M 

Owing to the Mobius invariance, the three-gluon state | m, M) , with the conformal 
weights m and M, can be written as a superposition of the Polyakov functions 

(\ -M-m+l 
P30_ \ 
Part/ 

with various integration contours L. By integrating the terms in M3 with derivatives of 
I po'j M) by parts and using the relations 

P30 \ . , w x-m+M / P30 \ 

— =-{M + m-l){p &0 ) — 

Pm'J KPso'J 

1 / \ -M-m+l 

1 „2 „ ofilT i\\ fl /„ . \-m+M-l / P30 \ 



( 1^2M ^ P ' ' " 2 9 °' P30 ' " 2(M ~~ 1} ) 93 (Po '° 



\P30'/ 

"(M + m - 1) (»"-^ M+1) („„.„)-— (M) """^ , 
one can obtain the recurrence relation for the function f m ,M = /m,M(pi, P2, P3', Po)'- 

M 13 ~ M 2s) fm,M = 1 _ 2M (^(M - 1) f m ,M+l + " ^ _ ~ "/m,M-l 

Due to its Mobius covariance f m ,M(pi, P2, P3] Po) can be presented in the form 

(\ m/3 
P12P23P31 \ r / \ 
2 2~ Jm,M\ X ) 1 
Pl0P20P30/ 

where x is the anharmonic ratio 

P12 P30 
X = . 

PlO P32 

By introducing the new integration variable 

' _ P12P30' 



PlO'P32 

we obtain the following expression for f m ,M( x ) '■ 

fmAx) = (X(l ~ X)) 2m/3 J dx\l -X)- M {X -X)-™^' 1 

This function satisfies the following differential equation 

M 2 12 {x) f mM (x) = M(M - 1) f mjM (x) 
and the recurrence relation 



x \ -M-m+l 



X 
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(MUx) - MUx)) f mM (x) = 
M(M + m-l) ( (m-M)(m-M + l) , \ 

1 _ 2M I ( M - !) fm,M+l( X ) + W^~l fm,M-l{x) I . 

The pair Casimir operators in the ^-representation are given below 

M i^) = ( J (! - 2 *) + - ^) ^ (f C 1 + *) + - » 

= - 2x) + x(l - ^ - 2) + x(l - x)d) , 

MUx) = -^^y (j(x - 2) + x(l - (|(1 + x) + x(l - x)d) 
and satisfy the relation: 

+ M^x) + M| 3 (x) = m(m - 1) . 

The function f m ,M(x) can be expressed for two different choices of the integration contour 
L through the hypergeometric functions: 

^■" {X) = r(M)r(2 - 2M) {X{1 ~ X))2 "" 3 xl ~ M ~ mF <~ m + 1 - M, 1 - M; 2 - 2M; x) , 



= r(2M)r7i + + !»-M) - I))2 '"" lM_mF(m + M - M; 2Mi *> • 

where 



_,, . . ^ ab x a(a + l)b(b+l)x 2 
F M ;e; *) = ! + -- + > - + .... 

Moreover, the above eigenvalue equation for f m ,M(x) is equivalent to the hypergeometric 
equation for F(a, 6; c; x): 

d 2 d 
x(l - x)—F + {c-(a + b + l)x) —F - abF = , 
dx l dx 

because 

(x(l - x)y 2m/3 x- M+m M 2 12 {x) (x(l - x)f m/3 x M ~ m = M(M - 1) + 

( d 2 d \ 

x \x(l- x)— + (2M - (m + 2M + l)x) — - M(M + m) . 
\ aar ax / 

In an analogous way, using the relation 

(x(l - x)y 2m/3 x~ M+m (M 2 3 (x) - M 2 3 (x)) (x(l - x)) 2m/3 x M - m = 
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-(M(2 - x) + (1 - m)x + (2 - (M(l - x) - m + x(l - x)d) 

Ob 

and extracting from it the hypergeometric differential operator with the coefficient chosen 
to cancel the term with the second derivative, we obtain that the above recurrence relation 
for fm,M(x) is equivalent to the following relation for the hypergeometric function Fm(x) = 
F(M,M ■ nr. 2 .1/ :.r): 

(m(m - 1) + M(2m — 1 — M) + 2(1 — m)(l — x)d + 2M (^~ m ) ^ Fjvf ^ 

. / (M — 1)(M + m)(M + m — 1) „ 2M „ .A 

= < M " "*> ( 2(2M-1)(2M + 1) X + ~ F "-^ j ■ 

One can verify its validity in the k-th order of the expansion in x by taking into account the 
algebraic identity: 

„,,, „,, rt 1 / M(M-m)\ (M + k)(M + m + k) 
(m + 2k)(m-l) + M(2m-l-M) + 2h-m+ ^ - - M * i 

(M — m)(M + m — 1) / (M — M(2M + fc-l)\ 

~ 2M - 1 \ 2M + fc + fc + 1 / ' 

Let us now return to the problem of finding the recurrence relations for the matrix V^(m) 
performing the unitary transformation between the M- and A-representations. As earlier, it 
is convenient to work with the function <p m (x) defined by the relation 

f m {x) = (x(\-x)y m '^ m (x). 

It satisfies the equation 

a x ^ m a m if m {x) = X(p m (x) , 

where 

a m = x(l-x)(id) m+1 . 
One can search the solution of the above eigenvalue equation as the linear combination 

M 

of the eigenfunctions of the operator M\ 2 . The coefficients Cjw(m) would coincide with the 
matrix V^(m), if the functions F^}(x) would be normalized. However, in accordance with 
above formulas we define F^{x) by the following expression: 

F£(x) = r(2M r[i t 1 +M)" M) {X(1 ~ X)r/3 / - M(X) = XM(1 " Xr F(m + M ' M; 2M; X) ■ 

The transition to the renormalized functions can be easily done. 

According to the above presentation the quantities F%}(x) satisfy the following equality: 

a^ m a m FZ (x) = ^M(M - 1)(M - m) ( ^^M-^M+l) K+l{x) " 2F ^~^) " 
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Because the right-hand side of this equality is zero for M — 0, 1 and m, one can restrict 
the summation over M in the eigenfunction tp m (x) to two series: M = r and M = m + r 
with r = 0, 1,2, .... However, in the first case the coefficient in front of F™(x) can not be 
calculated through the coefficient in front of Fq 1 (x). Due to this degeneracy of the equation 
one should introduce for r > 1 the more complicated function $™(x): 

The recurrence relation for these functions can be obtained by differentiating the relation 
for F^}(x) . In particular, for r = 1 we obtain 

a x . rn a m $r (x) = - m) ^ 1+ ^ m F ™ {x) _ 2F m (x) ) j . 

Thus, in accordance with the small- a; behaviour of the eigenfunctions ip(x), discussed 
above, we write the linearly independent solutions in the form: 



^HM) = E A r(A)iT(*)- 
k=i 

oo 

^ m) (^ A) = £ «(A) F?(x) + A™(A)$™ (x)) , 

k=0 

oo 

y4 m ^,A) = £ 7 r(A)Fr +m (x), 

k=0 

The coefficients a™, 7™ and A™ satisfy the recurrence relations: 

<AaHA)= (< + i(A)+/3r + i(A)^) k(* + l)(fc - ro+ 1) 

_i ( Qr . l(A) + ff _ i(A) a) (t _ _ 2)( , _ m _ i ^-wy ■ 

+A; (A; + l)(A;-m+l)A™ 1 (A), A™(A) = 1 ; 
1/, -x/, x (fc + 2m — l)(/c + 2m — 2) „ 

!A7r(A) = _ ( t + m - + ™ - 2 )( t - il^ri^t.W 

+ (* + m)(k + m + 1)(* + l) 7 r+l(A) • T„ m (A) = 1 . 

If we compare these relations with the derived above analogous recurrence relations for 
the coefficients of the expansion of Lp^ m \x) in the series over x, it is obvious that the factors 
in front of the corresponding quantities with the index k + 1 coincide. Furthermore, the 
quantities 7& and A^ are absent in the right-hand side of these relations contrary to the 
previous case, where the similar factors in front of , c& and dk are non-zero. 



31 



References 



1. L.N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 642; 

V.S. Fadin, E.A. Kuraev and L.N. Lipatov, Phys. Lett. B60 (1975) 50; 
Ya.Ya. Balitsky and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822; 
L.N. Lipatov, Sov. Phys. JETP 63 (1986) 904; 

2. V.S. Fadin, L.N. Lipatov, Phys. Lett. B429 (1998) 127; 

3. J.Bartels, Nucl. Phys. B175 (1980) 365 ; 

J. Kwiecinski and M. Prascalowicz, Phys. Lett. B94 (1980) 413; 

4. L.N. Lipatov, Phys. Lett. B251 (1990) 284; B309 (1993) 394 

5. L.N. Lipatov, |hep-th/9311037] Padua preprint DFPD/93/TH/70, unpublished; 

6. R.J. Baxter, Exactly Solved Models in Statistical Mechanics, (Academic Press, 
New York, 1982); 

V.O. Tarasov, L.A. Takhtajan and L.D. Faddeev, Theor. Math. Phys. 57 (1983) 163; 

7. L.N. Lipatov, Sov. Phys. JETP Lett. 59 (1994) 571; 

L.D. Faddeev and G.P Korchemsky, Phys. Lett. B342 (1995) 311; 

8. C. Montonen and D. Olive, Phys. Lett. 72 (1977) 117; 
N. Seiberg and E. Witten, Nucl. Phys. 426 (1994) 19; 

9. L.N. Lipatov, Recent Advances in Hadronic Physics, Proceedings of the Blois 
conference (World Scientific, Singapore, 1997); 

R.Janik and J. Wosiek, |hep-ph/9802TOO| , Crakow preprint TPJU-2/98; 

M.A. Braun, |hep-ph/9801352| , St. Petersburg University preprints; 

M.A. Braun, P. Gauron and B. Nicolescu, preprint LPTPE/UP6/10/July 98; 

M. Praszalowicz and A. Rostworowski, hep-ph/9805245| , Crakow preprint TPJU-8/98; 

10. R. Janik and J. Wosiek, Phys. Rev. Lett. 79 (1997) 2935; 

11. J. Ellis and N.E. Mavromatos, preprint OUTP-98-51P, |hep-ph/9807451 ; 



12. L.N. Lipatov, Perspectives in Hadronics Physics, Proceedings of the ICTP 
conference (World Scientific, Singapore, 1997). 

13. J. Maldacena, Adv. Theor. Math. Phys. 2 : 231 (1998), |hep-th/971l2T?D 



32 



